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Grain boundaries are commonly observed in carbon nanostructures, but their influence on thermal 
and electric properties are still not completely understood. Using a combined approach of density 
functional tight-binding theory and non-equilibrium Green functions we investigate electron and 
phonon transport in carbon based systems. In this work, quantum transport and thermoelectric 
properties are summarized for graphene sheets, graphene nanoribbons and carbon nanotubes with 
a variety of grain boundary types in a wide temperature range. Motivated by previous findings 
that disorder scatters phonons more effectively than electrons, a significant improvement in the 
thermoelectric performance for polycrystalline systems is expected. As the effect is marginally 
sensitive to the grain boundary type, we demonstrate that grain boundaries are a viable tool to 
greatly enhance the figure of merit, paving the way for the design of new thermoelectric materials. 


I. INTRODUCTION 

Graphene, the building block for various carbon nanos¬ 
tructures of all other dimensionalities, is the subject of 
exceptional scientific interest in the recent years mm- 
Today, graphene sheets, graphene nanoribbons and car¬ 
bon nanotubes are considered as promising candidates 
in the vast field of electronic, thermal and thermoelec¬ 
tric applications. However, macroscopic samples of those 
structures will most certainly appear in polycrystalline 
form, an issue which is neglected in most studies. Poly¬ 
crystallinity is an inherent structural impurity and very 
challenging to avoid. It is caused by growth kinetics or 
substrate imperfections, resulting in different graphene 
domains with a variety of crystallographic orientations 
[3]. The common domain interface, the grain boundary 
(GB), can be regarded as an one-dimensional array of dis¬ 
locations. In an infinite two-dimensional graphene sheet, 
this dislocation array can be approximated as a linear pe¬ 
riodic array [4] , whereas its periodicity or Burgers vector 
depends on the lattice mismatch. Besides unintentional 
polycrystallinity of the system as an inherent property, 
structuring of well-defined line defects has been demon¬ 
strated in graphene and opens new possibilities [5]. The 
resulting structures have been addressed both from ex¬ 
perimental and theory groups and can be classified by 
either their relative orientation angle or lattice mismatch 
[3 [SHE]. Despite a large number of experimental observa¬ 
tions, the implications on transport properties are not yet 
completely understood. It has been shown that heat flow 
in polycrystalline graphene depends on the specific GB 
structure [9] by effectively scattering phonons at the in¬ 
terface. These findings are encouraging to investigate the 
issue of thermoelectrics, which complements the recent 
studies focussing on either electronic or thermal trans¬ 
port mm- Additionally, we widen the scope of interest 
beyond two-dimensional graphene by including quasi-ID 
carbon allotropes, like graphene nanoribbons (GNR) and 
carbon nanotubes (GNT) in our studies. In fact, we show 
that GBs can not only tune a transport gap controlling 


charge currents [12] nn, but also significantly increase 
the thermoelectric figure of merit. Combined with the 
exceptional charge transport in carbon nanostructures, 
GBs seem to evolve as promising candidates in scatter¬ 
ing phonons to suppress the high thermal conductivity, 
which was found to reach up to 5000 Wm“^K“^ for free¬ 
standing single-layer graphene [14]. 

This work presents a summary of electronic and ther¬ 
mal transport properties in one-dimensional and two- 
dimensional carbon nanostructures with GBs in order to 
study the use of polycrystalline structures in thermoelec¬ 
tric materials. The paper is structured as follows. After 
a short introduction to motivate the topic in Sec. I, we 
will give an overview of the various structures and the 
theoretical framework for charge and thermal transport 
calculation in the next section of this paper. Sec. II. Re¬ 
sults will be discussed in Sec. Ill for each of the systems 
separately, i.e., graphene, GNRs and CNTs, and we will 
summarize and conclude subsequently. 

II. MODEL AND THEORY 

The polycrystalline nanostructures have been con¬ 
structed by joining two subsystems of different chiralities 
or crystal orientations. Based on those initial geometries, 
the minimum energy configurations of the structures were 
then determined using molecular dynamics simulation 
followed by density functional geometry optimization un¬ 
til the force convergence criterium of 0.02 eVA“^ was 
met. Hydrogen saturation of eventually di-coordinated 
carbons was neglected. In our analysis, the variety of 
GBs in graphene is limited by the need for a reasonable 
periodicity along the interface, otherwise the system size 
gets too large for ab initio calculations. This limitation 
is not apparent in the quasi-ID systems of GNRs and 
GNTs, but the choice of diverse chiralities is restricted 
by roughly matching the GNT diameters. 

For the analysis of quantum transport characteris¬ 
tics, the Green function formalism has been applied in 
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combination with the efficient density functional tight- 
binding approach (DFTB), as implemented in DFTB+ 
[isi [ig. This methodology allows the calculation of 
electron and phonon transport properties in the bal¬ 
listic transport regime based on the same theoretical 
footing ini HI]. Transmission spectras are obtained 
in the first place, which enable the calculation of elec¬ 
tron conductance, current-volt age characteristics, ther¬ 
mopower, thermal conductance and thermoelectric figure 
of merit. The here applied ballistic transport model ne¬ 
glects phonon-phonon and electron-phonon interactions, 
but the high intrinsic mean free paths for electrons and 
phonons at room temperature in carbon systems [8] val¬ 
idate this approximation. Furthermore, the model is 
superior to classical molecular dynamics approaches for 
thermal transport, as they lack quantum features like 
Bose-Einstein statistics. The high Debye temperature 
of about 2100 K in graphene systems [8] necessitates 
quantum calculations. The electron transmission func¬ 
tion Tq\{E) was obtained by the standard single-particle 
Green function formalism m- 


T,i{E) = Tr 


f£'effect 


G{E) 


ES-H 


-1 


( 1 ) 

( 2 ) 


The Hamiltonian H of the scattering region has been cal¬ 
culated using DFTB and semi-infinite unperturbed leads 
are assumed and incorporated by self-energy terms us¬ 
ing the decimation technique by Lopez Sancho et al 
m- Phonon transmission spectra Tph(ci;) were calcu¬ 
lated analogously based on the atomistic Green function 
method [18], 


rph(w) = Tr 




Qico) = 


luH-D 


(3) 

(4) 
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to the Landauer formula a = the temperature de¬ 

pendent Seebeck coefficient or thermopower 


S{T) 


Tli 

eTLo' 


(6) 


and the current-voltage characteristics for a source-drain 
voltage V applied between the two contacts 


nv) = ^J T,i(E) [ME,T) - ME + eV,T)] dE. 

(7) 

The thermal conductance consisting of an electronic 
contribution and the thermal conductance of the lat¬ 
tice /^ph, can then obtained by the phonon transmission 
function Tph 


^T t2/t \ ^ f rr ^ ^d/B(w,^) 

^ (7^2 - Ti/Lo) + ^ / - — -dw, 

l^el /Cph 

(8) 

where /b(<^, T) is the Bose-Einstein distribution function. 
In general, heat fiow is generated by lattice vibrations 
and free conduction electrons. However, in carbon nanos¬ 
tructures the electronic thermal transport is very limited 
[21] and accounts only for a few per mille of the total 
conductance. 

Einally, we compute the thermoelectric figure of merit 


ZT = 


aS‘^T 

K 


1 /^el 

(LoLa/q) - 1 ’ 


(9) 


which comprises of all previous quantities and charac¬ 
terizes the efficiency of the thermoelectric effect in the 
system. 


III. DISCUSSION OF RESULTS 


using the dynamical matrix D obtained from the mass- 
weighted force constant matrix, which can be evaluated 
from the second derivatives of the lattice potential en¬ 
ergy with respect to spatial displacements K = {kij} = 
d^UIdxidxj. This Hessian matrix has also been obtained 
by DFTB calculations and, in analogy of the electronic 
counter part, the leads account for self-energy terms cal¬ 
culated with the decimation technique. For further de¬ 
tails on the Green function formalism for electrons and 
phonons we refer to the corresponding literature or pre¬ 
vious publications mm- 

For a shorthand notation for several properties, we in¬ 
troduce the Onsager coefficients 

Ln{T) = I (5) 

with n G Mo, Fermi energy Ep and Fermi-Dirac distribu¬ 
tion /f(E,T). Based on the electron transmission func¬ 
tion Tel we calculate the electric conductance according 


A. Graphene 

The GBs in graphene sheets can be approximated as 
linear periodic arrays of dislocations [4]. We are using 
a GB classification in graphene proposed by Yazyev et 
al. [6] , separating into two classes corresponding to their 
matching vectors and (nR,mR). If exactly one 

matching vector fulfills the criterion {n—m) = {q £ Z), 
the GB is of class-H type. Otherwise it belongs to class- 
1. Due to a misalignment of allowed momentum-energy 
manifolds, class-H type boundaries introduce a transport 
gap, which can be approximated by Eq = hv-p/Sd = 
^ solely depending on the periodicity d [6]. The 
distinct behavior of both classes can be explained by 
transverse momentum conservation at the interface and 
an effective rotation of the Brillouin zone for the charge 
carriers passing the interface. We concentrate on two 
examples, one symmetric (class-I) and one asymmetric 
(class-H) GB, see Fig. For reference, the properties 
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of an unperturbed graphene sheet are calculated. As ex¬ 
pected, the asymmetric class-II GB exhibits an energy 
gap of about 1 eV, whereas the transmission spectrum of 
the symmetric class-I GB in Fig. [^A(i) is very similar to 
the pristine sheet. Phonon transmission shows a weak de¬ 
pendence on the GB type, but a slightly stronger phonon 
scattering can be identified for class-II GB with differ¬ 
ences most pronounced at very low (< 200 cm“^) and 
high phonon energies (> 1000 cm“^), see Fig. [^A(ii). 
This can be attributed to a stronger lattice deformation 
and buckling along the interface. The introduced asym¬ 
metry of electron and hole transmission leads to a sepa¬ 
ration of charge carriers, improving the thermopower S 
in those systems over ideal graphene. Electron holes ac¬ 
count for the Seebeck effect in class-II GBs, and therefore 
S' > 0, and give rise to an improved thermopower com¬ 
pared to the electron-dominated symmetric GBs. For 
the thermoelectric figure of merit ZT an enhancement 
by over three orders of magnitude in Fig. A(iv) is 
expected above room temperature for both types. In¬ 
terestingly, both GB types perform equally well, as the 
lower Seebeck coefficient S ~ 0.05 mVK“^ of the class-I 
GB is one order of magnitude smaller than for class-II 
type interfaces, S ^ 0.4 mVK“^ in Fig. [^A(iii), but 
gets compensated by an improved electron conductance. 


class-I, symmetric 



class-II, asymmetric 



(a) Examples of graphene GBs. The symmetric structure 
(2,1)1(2,1) corresponds to a class-I GB, whereas the asymmetric 
GB (5, 3) I (7,0) is class II. The corresponding matching vectors 
(nL,mL), are shown in black. 


GNRl, © = 24.0° GNR2, 0 = 16.0° GNR3, 0 = 11.0° 



(b) Examples of GBs in GNR structures, obtained by joining 
two '^1.7 nm wide GNRs with different inclination angles 0. 


B. Graphene nanoribbons 

As the gapless nature of graphene has negative ef¬ 
fects on its thermoelectric properties, spatial confine¬ 
ment leads to an intrinsic band gap which depends on 
size and edge structure. By parallel cutting polycrys- 
talline graphene sheets into ribbons, one expects struc¬ 
tures similar to those shown in Fig. i.e., at least two 
graphene nanoribbons of different crystallographic orien¬ 
tation with a GB at the interface. Gompared to graphene 
sheets, those structures omit periodicity parallel to the 
interface and reveal hydrogen termination at the edges to 
saturate dangling bonds and to remove particular edge 
state effects. The interface in between constitutes of an 
array of dislocations, either pentagons, heptagons, or oc¬ 
tagons. For one side of the system, an about 1.7 nm 
wide armchair GNR with 15 dimer lines in width was 
chosen. The ribbon across is rotated respectively by an 
angle 0, ranging from 11°to 24°. The angles are cho¬ 
sen in a way that the periodic supercell of the lead does 
not get unnecessarily large. By cutting the ribbon to the 
same width as the right part one produces irregular but 
periodic edge geometries. As a reference, the results for 
the unperturbed armchair ribbon are shown and, like all 
armchair terminated ribbons, it is semiconducting with 
a band gap of about 0.5 eV. As shown in Fig. 
the suppression of electron transmission is strongest in 
GNRl, which may be explained by the non-interrupted 
dislocation array, compared to the other configurations 
GNR2 and GNRS, where one hexagon in the interface is 
preserved. One also notices the electron-hole symmetry 


GNTI GNT2 GNT3 GNT4 



(c) Examples of GBs in GNT structures, obtained by joining two 
nanotubes with a similar diameter of about 0.55 nm. The 
chirality indices are given for each pair of GNTs. 


Figure 1. (Color online) Examples of polycrystalline car¬ 
bon nanostructures used for transport calculation. The GBs, 
obtained by DFT geometry optimization, exhibit pentagon 
(blue), heptagon (red), octagon (yellow) and nonagon defects 
(green). 


breaking, with an improved transmission for low energy 
electrons compared to respective holes. Interestingly, the 
effect of different GBs are not apparent in the thermal 
conductance. As expected, phonons get scattered at the 
interface, effectively lowering the transmission by roughly 
50%, but the spectrum is mainly independent on the ex¬ 
act geometry, see Eig. [^B(ii). This results in about one 
tenth of the thermal conductance of an ideal ribbon. Ef- 
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GNRl — GNR2 — GNR3 


(G) GNT 


— GNTl — GNT2 — GNT3 — GNT4 


0.001 



T(K) 


T(K) 


^Oo 

T (K) ^ 


Figure 2. (Color online) Calculated electron transmission (i) and phonon transmission (ii) at room temperature, as well as 
the Seebeck coefficient (iii) and thermoelectric figure of merit (iv) as a function of temperature for graphene (A), graphene 
nanoribbons (B) and carbon nanotubes (C) with different GB types, see Fig. Black curves represent the monocrystalline 
system for reference. 


fects on the electronic and thermal properties combined 
give rise to a significantly improved thermoelectric fig¬ 
ure of merit, shown in Fig. [^B(iv). At room temper¬ 
ature, we calculated an increase by three to four orders 
of magnitude, topping at around ZT 0.1 for GNR3 at 
700 K. The zero in the ZT graph of GNR3 can be associ¬ 
ated to a change of the majority charge carrier type from 
holes to electrons with increasing temperature, implying 
a sign change of the Seebeck coefficient S at 210 K. In 
fact, S is negative with values around S ~ —0.4 mVK“^ 


for all configurations, in contrast to S' > 0 in the ideal 
GNR. Those findings are particularly promising as all 
three samples show the same qualitative performance, 
suggesting a robust and efficient method for enhancing 
the thermoelectric effect. 
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C. Carbon nanotubes 

Next to the planar carbon structures, their rolled-up 
form is equally promising. Here we will shed light on 
the electric and thermal properties of carbon nanotube 
heterojunctions. Carbon nanotubes of different chirality 
but similar diameter can form a junction by exhibiting 
a couple of dislocation defects, see Fig. Such struc¬ 
tures can evolve if two individual tubes eventually grow 
together, or by a change in the growth parameters 
[25] . First interest in application of intra molecular nan¬ 
otube junctions is the use as hetero junction diodes [26] . 
We studied four different CNT heterojunctions, each of 
them with a diameter of about 0.55 nm. In all samples 
one part is made of a (7,0)-CNT, which is complemented 
by a second CNT with aligned tube axis, see Fig. [T^ 
Kinked carbon nanotube junctions are not considered in 
this work. Notice, that the structure CNT3 in Fig. [T^ 
shows no chirality change, but a twist around the tube 
axis. The twisting angle is given by (or multiples of) 
0 = ^ 12.8° due to the seven-fold symmetry. This 

twist deforms the lattice and increases stress at the inter¬ 
face, but does not lead to defect formation or reconstruc¬ 
tion. In the remaining three structures CNTl, CNT2 and 
CNT4, the chirality indices (5,3), (4,4) and (6,1) have 
been chosen, all yielding a diameter close to its coun¬ 
terpart of 0.55 nm. Primarily one observes pentagon- 
heptagon pairs along the circumference, and one octagon 
in CNT2. The properties of the ideal and semiconduct¬ 
ing (7,0)-CNT are calculated as a reference. Fig. |C(i) 
shows the energy gap in the electron transmission ap¬ 
pearing in all of the different heteroj unctions. Most no¬ 
tably, the twisted carbon nanotube, CNT3, shows only 
little change in the electron transmission where the slight 
lattice perturbations smear out the perfect transmission 
steps of an ideal nanotube. Electrons in CNTl and CNT4 
are suppressed equally and transmission is roughly cut in 
half. The heteroj unction in CNT2 shows the strongest 
scattering of electrons, effectively lowering electron con¬ 
ductance. We attribute this to the strong deformation 
caused by the octagon defect. A similar picture can be 
drawn for phonon transport, see Fig. |^C(ii). As no de¬ 
fects are formed in CNT3, lattice vibrations are weakly 
affected and scattering happens only for high frequency 
phonons. But for heteroj unctions of different chiralities, 
phonon transmission gets significantly suppressed at all 
frequencies. Interestingly, CNT4 with only one defect 
pair is superior to structure CNTl in blocking phonons 
and the defect configuration in CNT2 scatters strongest. 
Consequently, thermal conductance is decreased up to a 
factor of five for the CNT2 structures compared to the 
ideal nanotube. Quantitatively, this is very similar to 
the case of GNRs discussed above. In terms of figure of 
merit, the thermoelectric properties can be significantly 
improved but strongly depend on the exact structure, see 
Fig. [^C(iv). At room temperature, a gain in calculated 
ZT as high as two orders of magnitude for the struc¬ 


tures CNTl and CNT2 can be reported. Note that the 
twisted structure without chirality change, CNT3, shows 
a decreased ZT by a factor of ten. This is due to the 
loss of electron-hole asymmetry as the slight distortions 
smear out the sharp band edges. Thereby the Seebeck 
effect is strongly reduced. With the right choice of CNT 
chiralities, the Seebeck effect can be increased to values 
around S ~ —0.6 mVK“^ for a broad temperature range 
and a substantial enhancement of ZT is expected. 


IV. CONCLUSIONS 

A material with high thermoelectric efficiency has 
to show properties of an electron crystal and phonon 
glass at the same time, i.e., blocking thermal transport 
while maintaining electric current. Our results show 
that lattice imperfections forming at the interface of 
two carbon systems with different crystallographic ori¬ 
entations scatter heat stronger than charge carriers. In 
studying the effect of various GB types on the electron 
and phonon transport properties of polycrystalline car¬ 
bon nanostructures, we predict improved thermoelectric 
properties. The odd-membered rings at the boundary 
break the bipartite symmetry of the lattice resulting in 
an electron-hole asymmetry in the electron transmission 
spectrum. This charge carrier separation is advantageous 
for the thermoelectric effect, increasing the thermopower. 
This, combined with suppressed thermal transport due 
to phonon scattering at atomic dislocations, in particu¬ 
lar for mid to high energy phonons, generates a substan¬ 
tial improvement in the figure of merit. Our calculations 
show for nearly all tested configurations an enhancement 
of several orders of magnitude in ZT at room tempera¬ 
ture. The low sensitivity of this effect on the GB type 
heavily facilitates any experimental realization. We also 
note that there is even more room for improvement of 
ZT by shifting the chemical potential. This opens up 
promising prospects for the use of polycrystalline carbon 
nanostructures in thermoelectric applications. 
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